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ABSTRACT 

We present an analysis of the AGN broad-line regions of 6 powerful radio galaxies at z~2 (HzRGs), which is part of a study of a sample 
of 50 HzRGs with rest-frame optical imaging spectroscopy obtained at the VLT. In 6 galaxies we detect luminous (»S?(Ha)=fewx 10 44 
erg s~ 1 ), spatially unresolved, broad (FWHM> 1 0,000 km s~ 1 ) Ha line emission from the nucleus (HaBLRs), consistent with broad- 
line regions of supermassive black holes with masses of fewx 1O 9 M0 and accretion luminosities of a few percent of the Eddington 
luminosity. In two galaxies we also detect HjSBLRs, suggesting relatively low extinction of Ay ~ 1 mag, which agrees with constraints 
from X-ray observations. Overall, we find HaBLRs in ~20% of the galaxies where Ha is observed. By relating black hole and bulge 
mass, we find a possible offset towards higher black-hole masses of at most ~ 0.6 dex relative to nearby galaxies at a given host 
mass, although each individual galaxy is within the scatter of the local relationship. If not entirely from systematic effects, this would 
then suggest that the masses of the host galaxies have increased by at most a factor ~ 4 since z~ 2 relative to the black-hole masses, 
perhaps through accretion of satellite galaxies or because of a time lag between star formation in the host galaxy and AGN fueling. We 
also compare the radiative and mechanical energy output (from jets) of our targets with predictions of recent models of "synthesis" or 
"grand unified" AGN feedback, which postulate that AGN with similar radiative and mechanical energy output rates to those found 
in our HzRGs may be nearing the end of their period of active growth. We discuss evidence that they may reach this stage at the same 
time as their host galaxies. 
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1. Introduction 

One of the most remarkable discoveries in extragalactic as- 
tronomy in the past decade is the close relationship between 
the masses of supermassive black holes and several properties 



of their bulges, such as stellar mass (Magorrian et al. 1998 
Haring & Rix 20 04), ste llar velocity dispersion ("M-a rela 



etal.|2 002; 
•aham et al. 



tion" Gebhardt et al. 2000; Ferrarese 2002; Tremaine 
Ferrarese 2002|, and the concentration parameter ( jGraham 



2001). Theoretical arguments suggest that the most fundamental 
of these may be a relationship that with the binding energy of 
the bulge (Hopkins et al. 2007; Young eT et al.|20 08 ) or perhaps 
with the dark-matter halo ( |Booth & Schaye|2010| l. 

A popular way to interpret these relationships is by pos- 
tulating that supermassive black holes may regulate their own 
growth, as well as that of their host galaxies, by coupling a small 
fraction of their energy output to the surrounding gas. This heats 
and unbinds significant fractions of the gas, by inhibiting star 
formation ("AGN feedback" e.g. pF& Rees||1998[ |Springel | 
et al. |2005 HBower et al.|20061|Croton et al.|2006l|Cattaneo et al.| 
2007] I Ciotti & Ostriker l |2007[ JSijacki et al.||2007| [Merloni &| 



Heinz|2008[|Ciotti et al. 



2009; Fanidakis et al.|20091l. The scal- 



* Based on observations collected at the Very Large Telescope of 
ESO. Program IDs 070.A-0545, 070.A-0229, 076.A-0684, 079.A-0617, 
081A-0468,381A-0541, 082.A-0825, 083.A-0445, 085.A-0897 



ing relations of black hole mass with the bulge parameters are 
a major observational constraint to gauging the "efficiency" of 
this process, the fraction of the energy output of the AGN that 
is being injected into the interstellar medium of the host galaxy 
(typically a fewx 0.1%, e.g. [DTMatteo et al.||2005| |Merloni & 



|Heinz|2008[|Shankar et al.|2009[|Kording et al.|2008| l 

A large number of observational studies have been dedicated 
to finding direct observational constraints for the (co-)evolution 
of black holes and their host galaxies across cosmic time (e.g. 



Shields et al. 2003; Walter et al. 2004; McLure et al. 2006; Peng 


et al.||2006 IShields et al. 2006[ |Woo et al. 


2006[ ITreu et al. 


2007[ |Alexander et al. 2008 Coppin et al. 


2008 Shen et al. 



redshifts z>2 are particularly interesting, since these redshifts 
correspond to the epoch when powerful AGN were most abun- 
dant (Pei 1995) and when massive galaxies evolved most quickly 
(e.g., Caputi et al. 2006). However, since we know little about 
the underlying physics of how AGN interact with their environ- 
ments and whether feedback occurs in AGN and host galaxies of 
all types at these redshifts, it is difficult to link black-hole prop- 
erties and feedback directly. This is particularly true in light of 
recent galaxy evolution models, which postulate that supermas- 
sive black holes may interact with their environments at different 
stages of their evolution in different ways, where thermal feed- 
back may dominate at some epochs, and mechanical feedback 
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may at others (|Churazov et al.|2005||Croton et al. 


2006HSijacki 


et al.|2007 
et al |2009 


Merloni & Heinz 2008; Fanidakis et a 


.120091 ICiotti 


K 



Using integral-field spectroscopy of the ionized emission- 
line gas in high-redshift radio galaxies at z>2 we recently iden- 
tified energetic, kpc-scale outflows of ionized gas correspond- 
ing to significant fractions of their interstellar medium (a few 
xlO 10 M ; |Nesvadba et al.|2006l |2007[ |2008b} >. These ionized 



gas masses are comparable to the cold molecular gas masses in 
these galaxies (e.g., Papadopoulo s et al.|2000| |De Breuck et al.| 
|2003| [2005| [Ne svadba et al. 2009), suggesting that a significant 
fraction of the interstellar medium participates in the outflows. 
Outflow velocities are of order of the escape velocity of massive 
galaxies. These may be the observational signatures of energetic 
AGN-driven winds, making the black holes of HzRGs particu- 
larly interesting targets. 

Unfortunately, nuclear broad-line regions (BLRs) in HzRGs 
are very rarely observed, because the central regions are ob- 
scured in these type II AGN. This has significantly hampered 
observational efforts to characterize their black-hole properties. 
Most notably, at redshifts z>2, BLRs have only been anecdo- 
tally reported for 4 HzRGs (|Rottgering et al.|1997| |Larkin et al 
{Nesvadba et al.||2006| |Humphrey et al.||2007| >, and their 



2000 



black- 



lole properties have not yet been fully discussed. McLure 



|et al.| ( |2006| l studied the redshift evolution of the black-hole 
bulge mass relationship for powerful radio-loud AGN out to z<2 
(with 6 radio galaxies at 1.5<z<2.0 and one quasar at z=2.0), 
but relied on the BLRs of quasars for the black-hole properties 
and on radio galaxies for the host properties, assuming that the 
only difference between radio-loud quasars and radio galaxies in 
a given redshift range was in their orientation. This has not yet 
been proven for high-redshift AGN, where evolutionary effects 
may be at least as important (e.g. Hopkins et al.|2006 i. 

Here we discuss the properties of supermassive black holes 
in 6 powerful radio galaxies at z~2-2.5 (HzRGs), in which we 
detect broad Ha line emission with FWHM> 10,000 km s -1 
originating in a compact region near the nucleus. The similarity 
to broad-line quasars suggests this is nuclear broad-line emission 
from the AGN. The emission lines are luminous, but less so than 
in powerful, optically selected quasars, making it possible to 
study AGN and host properties for the same objects. We discuss 
black-hole masses and accretion rates, and use the stellar mass 
estimates of Seymour et al. (2007) to place these galaxies on the 
black-hole bulge mass relationship. We also compare radiative 
and mechanical energy output rates with predictions from recent 
AGN feedback models (Churazov et al.|2005|JSrJacki et al.|2007 



Merloni & Heinz 2008; Fanidakis et al. 2 009} , which suggest 
that black holes with similar properties as those in our targets 
may be in a critical transition phase from radiatively efficient to 
radiatively inefficient accretion. 

Throughout the paper we adopt a flat Ho= 70 km s~ 1 Mpc~ 1 
concordance cosmology with Q.m = 0.3 and Q.a =0.7. 



icant fraction of the known radio galaxies with confirmed red- 
shifts at z>2 observable from the VLT and at suitable redshifts 
for near-infrared observations. All are type-2 AGN, classified by 
narrow rest-frame UV emission lines and a high Hell/CIV line 
ratio. The results of the full survey will be discussed in a series 
of upcoming papers. Here we focus on the nuclear properties of 
the 6 galaxies that have broad Ha line emission with proper- 
ties consistent with originating in AGN broad-line regions. We 
defer a full analysis of the spatially-extended line emission to 
a subsequent paper, but emphasize that the present galaxies are 
consistent with having outflows very similar to those described 
by |Nesvadba et al.| ( |20^6l|2007]|2008b| ). 

Most data were obtained with the H+K grating giving a 
spectral resolving power R=1500 and a spectral coverage be- 
tween 1.45/im and 2.4/im. Total on-source exposure times were 
between 45 min and 6 hrs per source. MRC1 138-262 was ob- 
served in the H and the K band individually, with R^ 3000 and 
4000, respectively, and with on-source exposure times of 75 and 
140 min, respectively (The data of MRC 1138-262 have previ- 
ously been described by [Nesvadba et al.||2006| but without a 
dedicated analysis of its nuclear properties.) We reduced all data 
sets with our custom data reduction software, which is optimized 
to study high-redshift galaxies. For details of the data reduc- 
tion procedure, see, e.g., |Nesvadba et al.| ( |2006| |2008a[ ). Data 
were calibrated using nearby bright stars observed at similar air 
masses to the targets at the beginning or end of each observa- 
tion. We used the same stars to monitor the size of the seeing 
disk with SINFONI and the same setup as for the data. Typically, 
the seeing was between 0.5" and 1.5", corresponding to a spatial 
resolution of 4—12 kpc at a redshift z^2.5. 

Three of our sources have already been known to have broad 
nuclear Ha line emission, but their black-hole properties have 
not yet been discussed in the literature. |Larkin et al.| ( |200"0| > dis- 
covered a broad nuclear Ha line in MRC2025-218 at z=2.6. 
[Nesvadba et al.| ( |2006| l identified a broad Ha emission line in 
MRC1 138-262 at z=2.2. [Humphrey et al.| ^2007\ report on the 
discovery of broad nuclear Ha emission in MRC1558— 003 at 
z=2.5. For MRC1558-003 and MRC2025-218 our measure- 
ments are consistent with the previous results. For TXS1113- 
178, Rottgerin g et al.| ( 1997| > find a broad CIV component con- 
sistent with the FWHM of our Ha line measurement. lEales &l 
|Rawlings] ( |1996| l report the detection of a "broad" Ha line in 
MRC0156-252; however, their spectrum has low signal-to-noise 
and they give an FWHM estimate of about 4000 km s . This 
is consistent with the integrated extended, spatially-resolved 
line emission in our data, and significantly narrower than the 
FWHM=12500 km s -1 we find for the broad line. In the fol- 
lowing we use our new measurements, adopting the results pub- 
lished in |Nesvadba etaT[ ( |2006l l for MRC1 138-252. A few broad 
lines were previously detected in radio galaxies at lower red- 
shifts, z~ 1 -1.5 ( |Economou et al.|1995[|Willott et al.|20"00l >. 



2. Observations and sample selection 

All data were taken with the near-infrared integral-field spectro- 
graph SINFONI ( |Tecza et al.1|2"0"00"l |Bonnet et aLl[2004) > at the 
Very Large Telescope of ESO in the H and K bands. SINFONI 
is an image sheer giving a field of view of 8" x 8" with a pixel 
scale of 250 mas in the seeing-limited mode. These observations 
are part of an ongoing survey of HzRGs at z>2. Our sample 
currently consists of 50 powerful radio galaxies at z>2, which 
were selected according to radio power and radio size. All tar- 
gets have steep-spectrum radio sources. They include a signif- 



3. Results 

Figure [T] shows the integrated spectra of our targets centered on 
the position of the brightest continuum emission. Spectra were 
extracted from 2"x2" box apertures (corresponding to about 
twice the FWHM of the seeing disk). Broad emission-line com- 
ponents are apparent in all six galaxies, and are superimposed 
by luminous, narrow components of Ha and [NII]AA6548,6583. 
MRC 101 7-220 is at a redshift where the red wing of the broad 
line falls into a spectral region with very low atmospheric trans- 
mission between the H and the K band, so that we only detected 
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the blue wing of the broad line and the line center, but missed 
most of the red wing. 

Owing to the relatively broad widths of the narrow lines of 
typically FWHM^800 km s -1 , [Nil] and the narrow Ha com- 
ponent are blended in all galaxies. The kinematics of these lines 
is consistent with those in the extended emission-line regions, 
suggesting they are part of the same outflows as seen in projec- 
tion along the same line of sight as the AGN. We will discuss the 
extended line emission in a forthcoming paper. 

3. 1. Broad Ha line emission 

To derive emission-line properties, we removed the continuum 
emission and fitted all broad and narrow line components with a 
set of weighted Gaussian profiles. To distinguish between broad 
and narrow components, we attributed zero weights to wave- 
lengths where the narrow components (or night-sky lines) domi- 



natq^J We subtracted the broad components, and fitted the resid- 
uals with lines whose redshifts and line widths were tied to the 



values observed in [OIII]A5007. |Nesvadba et aL] ( |2008b| ) find 
that, in the extended emission-line regions of HzRGs, all bright, 
narrow emission lines typically have very uniform kinematic 
properties, at least to the precision of the observations, and we 
find the same in the present analysis. 

Broad Ha lines typically have fluxes Fn a ~ 10~ 14 — 10~ 15 
erg s" 1 cm" 2 and FWHMs- 10000- 14000 km s -1 . Values for 
individual targets are given in Table [T] We did not correct for 
Galactic foreground extinction, which is negligible at the ob- 
served near-infrared wavelengths. One broad and one narrow 
component were sufficient for all galaxies except MRC0156- 
252, for which we needed an additional Ha component of in- 
termediate width, FWHM=7480 km s _1 and an integrated line 
flux of 1.4 x 10~ 14 erg s _1 cm~ 2 (Figure]^}. 

All data sets are deep enough that we would have detected 
Ha lines with similar luminosity and width, which suggests that 
HzRGs with HaBLRs may have different AGN properties or 
may be more highly extincted than other HzRGs. To quantify 
this we measured the r.m.s. at the spectral position of Ha in each 
of the cubes where Ha falls into the band, but was not detected. 



1-19 



erg s cm , and up 
1 x 10~ 18 erg s _1 cm~ 2 for the shallowest data sets. We 



We find a typical r.m.s. of a fewx 10 
to 

used Monte Carlo simulations to estimate 3(7 limits on the in- 
tegrated flux of these lines, assuming two line widths which are 
representative of the measured widths. For 1000 iterations, we 
can exclude line fluxes brighter than 8.7 x 10~ 17 erg s _1 cm -2 
and 2.6 x 10~ 16 erg s _1 cm -2 at 3d significance for line widths 
of FWHM=9500 km s" 1 and FWHM=12500 km s -1 , respec- 
tively. This is about an order of magnitude fainter than the fluxes 
of the detected broad Ha lines. 



3.2. Broad Hp line emission and extinction 

In the two galaxies with the deepest observations we also de- 
tected a broad H/3 component in the H-band. Figure [3] shows 
the broad Hj3 line emission in MRC0156-252 and MRC1558- 
003. TXS1113-178 and MRC1138-262 have more shallow data 
which are not sufficient for probing Hj3 line emission in a simi- 
lar way. In MRC0156-252 and MRC1017-220, Hj3 falls outside 
the near-infrared atmospheric windows. 



1 For MRC1017-220 we only fitted the blue wing, so that we are not 
able to account for line asymmetries. By analogy with the other targets, 
and given the large astrophysical uncertainties, we do not think that this 
has a major impact on the results. 



Although the broad Hj3 lines are not intrinsically faint, 
they are difficult to detect due to their large widths com- 
pared to the resolution, imperfections in the night-sky line 
subtraction, superimposed narrow-line components, including 
[OIII]AA4959,5007, and the shape of the continuum. We there- 
fore did not fit the redshift and widths of these lines individually, 
but used the values found for Ha to identify the lines and to es- 
timate their fluxes. In Figure[3]we show the expected H/3 profile 
for both galaxies. By scaling this spectrum we derive Ha/H/3 
flux ratios of F(Ha)/F(Hj3)^4-5 for both galaxies. By compar- 
ing our spectra with the Fell quasar template of Iwamuro et al.| 



( |2002| l we can rule out strong Fell contamination at the wave- 
length of Hj3. 

For an intrinsic line ratio of Ha/Hj3 =3.1 and a galactic ex- 
tinction law, the observed Ha/Hj8 ratios correspond to an ex- 
tinction of Ay ~1 mag. [Dong et al. ([2008 ) report very similar 
Ha/H/3 ratios in radio-loud AGN at low redshift, Ha/H/3 =3.4. 

Equally low extinctions are also suggested by comparing the 
observed Ha and X-ray fluxes. If most of the Ha line emis- 
sion is due to photoionization of gas in the central few parsecs 
of the AGN, then simple energy conservation arguments imply 
that the Ha line emission cannot be more luminous than the X- 
ray emission. In fact, observations of nearby AGN suggest X- 
ray/Ha ratios of 3-20 ( Ward et al.| 1988 Imanishi & Terashima 
2004). X-ray l uminosities are availa ble in the literature for two 
of our sources. |Overzier et al.| ( |2005| l m easure Jzfy = 2 x 10 45 erg 
s 1 at 2-6 keV for MRC0 156-003, and |Carilli et"ai"l l |2002) find 
_gf z = 4 x 10 45 erg s" 1 at 2-10 keV for MRC 1138-262. These 
X-ray and Ha luminosities given in Table [2] imply a ratio of X- 
ray to Ha emission of ~ 6 for both sources. This is within the 
range of values suggested by Ward et aL] ( |1988| l and Imanishi 
|& Terashima| ( |2004| l, and is fully consistent with relatively low 
extinction along the line of sight into the AGN. | Alexander 
|et al.| ( |2008[ ) give similar arguments for submillimeter-selected 
quasars and point out that this also suggests that scattering is un- 
likely to play a large role for the observed Ha line fluxes. The 
same applies to our HzRGs. 

The same comparison for two HzRGs with similar X-ray and 
near-infrared data ( |Overzier et al.|2005) [Nesvadba et al.|20 08b) 
suggests that the HaBLRs of HzRGs may generally suffer much 
higher extinction. |Overzier et al. ([2005 ) measure _2x = 1 x 10 44 
erg s" 1 for MRC0406-244 a t z=2.4 and ££ x = 6 x 1 44 erg s" 1 
for TXS0828+193 at z=2.6. ( |Nesvadba et al.|2008b| ) did not re- 
veal any HaBLRs but allow us to place upper limits on their 
observed Ha luminosity of 2 x 10 erg s . This corresponds 
to much smaller fractions of their X-ray luminosity, 1/20 for 
MRC0406-244 and 1/100 for TXS0828+193. It may indicate 
that extinction depends strongly on the individual sightline into 
the AGN. This is plausible if the spatial distribution of the ob- 
scuring material is inhomogeneous and clumpy. 



4. The supermassive black holes of HzRGs 

4. 1. The frequency of HaBLRs in HzRGs 

Previous discoveries of HaBLRs in HzRGs have only been 
anecdotal ( Larkin et al.|2000| [Nesvadba et al.|20 06, Humphrey 



[et"aT1|2007 



see also Econom ouet al.||1995[|Willott et al.||2000 
at somewhat lower redshifts, z<1.5). Although we did not se- 
lect a complete sample in the usual statistical sense, our targets 
include a significant fraction of the known HzRGs at these red- 
shifts, and may therefore allow a rough assessment of how fre- 
quent HaBLRs may be among the most powerful radio galaxies 
known at z>2. Having 27 targets at redshifts where HaBLRs 



4 



Nesvadba et al.: The black holes of HzRGs 



fall into the atmospheric windows and adding 4 HzRGs at simi- 
lar redshifts of Humphrey et al. ( 2007| l which are not in our sam- 
ple, we find that HaBLRs were detected in 6/3 1 possible cases, 
corresponding to ~ 20% of our sample. 



4.2. The nature of HaBLRs 

HzRGs typically have luminous, extended line emission with 



|Jarvis & McLroe] ( |2002" 



comple x kinematics. For example, Nesvadba et al. (2006 2007 
|2008b| ) find line widths of up to FWHM~4000 km s^Tn the in- 
tegrated spectra of HzRGs, which are due to spatially extended 
gas with large velocity gradients and line widths, and could 
mimic a BLR (see the example of MRC0156-262,Sj2|. In our 
targets the Ha emission-line morphology shows a luminous, nu- 
clear, and compact Ha emission-line component in addition to 
the more diffuse extended (and narrower) line emission. These 
broad components have no equivalent in [OIII]A5007, unlike 
the emission from the diffuse gas, where the spectr al profiles 
of Ha and [OIII]A5007 are typically very similar (Nesvadba 
|et al.|[2007] [2008b ). The absence of broad components of for- 
bidden lines is among the prime characteristics of BLRs where 
the electron densities are high enough that forbidden levels are 
collisionally de-excited. 

We take this as evidence that the HaBLRs of our targets in- 
deed trace gas in the vicinity of the supermassive black hole, and 
adopt the common interpretation that this gas is approximately 
in virial equilibrium and is an adequate tracer of the black hole 
mass. Accretion rates are relatively low in our sources 4.4 1 
well below the values for which we may expect strong disk 
winds ( |Kurosawa & Proga|2009| |Proga & Kallman|2004| . 



4.3. Virial mass estimates 

A number of recipes have been proposed to estimate black-hole 
masses from the line widths and luminosities (e. g., |McLure & 
Jarvis|2002||Vestergaard|2002||Onken et al.|2004}|Peterson etaL 



2004| |Vestergaard & Peterson||2006[) , among which we choose 
the approach of |Greene & Ho| ( |2005| l which relies on the prop- 
erties of the Ha line only. This is the obvious method of choice, 
because Ha is the only line that is well detected in our targets. 

|Greene & Ho| ( |2005[ ) give a simple scaling for black-hole 
mass, Mbh, as a function of the FWHM and luminosity of the 
Ha line, M BH = 2 x 10 6 jSf£ 55 x FWHM\ 06 , where FWHM 3 
and Ha luminosity ^2 are given in units of 10 3 km s _1 and 



10 erg s . With this scaling and the luminosities and FWHMs 
given in Table [T] we find a relatively narrow range of black-hole 
masses of a few x 10 9 M and up to 2 x 10 10 M . Results for 
individual galaxies are given in Table [2] 

Applying a method calibrated for quasars to radio galaxies 
can potentially introduce important systematic uncertainties, in 
particular owing to inclination and extinction. Section 13. 2| sug- 



gests that ext inction does not likely play a major role. |Greene 
& Ho (2005 i find an average observed ratio of HafHf} «3.5, 
whereas we find Ha/Hj3 «4. Since black-hole mass depends 
only weakly on the Ha luminosity, this difference is negligible. 

Orientation is perhaps a bigger concern, if our radio galax- 
ies are on average more strongly inclined relative to the line of 
sight than quasars ( e.g., Antonucci|1993| l. Greene & Ho ( 2005 ) 
do not find any significant difference between radio-quiet and 
radio-loud quasars, but it is likely that their radio-loud sample 
consisted mostly of flat-spectrum quasars (although they do not 
comment on the radio spectral index). 



) compare quasars with flat and steep- 
spectrum radio sources and find on average a factor 1.3 broader 
lines in steep-spectrum quasars similar to our HzRGs, which 
would correspond to a bias towards higher masses by a factor 
1.8. Of course we cannot exclude that HzRGs may on average 
be even more strongly inclined than steep-spectrum radio-loud 
quasars. If, for example, the BLR in such a quasar had an incli- 
nation of ~ 45° relative to the line of sight, and that of a radio 
galaxy an inclination of ~ 80°, then we would expect to overes- 
timate the black hole mass by another factor ~ 2, which is not 
very different from other astrophysical uncertainties (see also, 
e.g., |Collinetal.|2006) . 



Orientations near the line of sight are also disfavored by rela- 
tively uniform and moderate fractions of the luminosity of the ra- 
dio core relative to the lobes of 0.7% to 10% ( |Carilli et al.|1997[ 
|Pentericci et al.|2000|l, which suggest orientations between 45° 
and 60° (e.g.. Carilli et al. 1997| l. This also suggests that orienta- 
tion probably does not introduce a major bias. Moreover, it is not 
clear whether the unified model is appropriate for high-redshift 
AGN, since evolutionary effects may at least play as large a role 
at high redshift (e.g. Hopkins et al. 2006 see also Sj6]). When 
necessary in the following we state both mass estimates, the one 
with and the one without inclination correction. 



4.4. Eddington luminosities, radiative and mechanical 
luminosities 

The accretion luminosity of a black hole is the second fundamen- 
tal parameter we estimate from our data, by relating the bolo- 
metric luminosity with the maximally possible accretion lumi- 
nosity for a spherically symmetric flow (Eddington luminosity). 
The Eddington luminosity, Jz?Edch of a SMBH sc ales with bla ck- 
hole mass as Jz? Edd = 1.3 x 10 38 M/M@ erg s" 1 ( |Rees|l984| . In 
Table [2] we list the Eddington luminosities corresponding to the 
black-hole masses estimated in 94.31 



The bolometric luminosity can be estimated from our data. 
|Greene & Ho| ( |2005| l calibrated the broad Ha luminosity as 
a function of the continuum luminosity at 5100 A, ^ua = 
5.25 x 10 42 ^510044. and we adopt -^boi = 9 A £? 5m to translate 
the continuum luminosity at 5100 A to the bolometric luminosity 
(Kas pi et al.|2000[ ). With the black-hole properties listed in Table 
|2J we find accretion luminosities of 0.02-0.07 ££ Edd (0.04-0. 14 
J^Edd f° r inclination-corrected black-hole masses, §4.3| l. This re- 
sult is not dominated by extinction. Ha/Hj3 line ratios are not 
much higher than those of quasars ( |Greene & Ho 20051, and 
Eddington luminosity and bolometric luminosity only weakly 
depend on ^£ua with powers of 0.55 and 0.86, respectively. 
Therefore extinction is unlikely to lower the Eddington ratios 
by more than ~ 10%. 



4.5. Kinetic luminosity of the radio source 

In Sj6] we rely on the radiative and mechanical energy output 
rates of the supermassive black holes. This requires constrain- 
ing the kinetic jet luminosity, which is inherently difficult, be- 
cause the observed radio luminosity is not a good tracer of the 
energy content of the relativistic jet plasma. Consequently, jet 
kinetic powers are only known at an order-of-magnitude level 
or worse (e.g., Bicknell 1995). We adopt two different methods 
to estimate the kinetic power in our sources. Both methods rely 
on low-frequency fluxes, which we estimate from extrapolating 
along the power law of the radio synchrotron emission based 
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on multi-frequency observations listed in the NED extragalactic 
data base. 

First, we use the calibration of Birz an et al.| (|2008 ), who es- 
timated the energy necessary to inflate cavities in the diffuse 
hot gas of massive galaxy clusters. This is a very direct em- 
pirical estimate of the jet kinetic power, but was derived for 
radio sources that are 2-3 orders of magnitude less powerful 
than those in our targets. Second, we rely on the results of 
|Willott et al.| (|T999), who propose a scaling of the kinetic power 
as a function of the monochromatic radio flux at 151 MHz, 

■&kin = 3 x 10 45 f^j 2 -&fsi28 er S s_1 - ^ ere -^151,28 indicates the 
measured flux density at 151 MHz in units of 10 28 W Hz 1 sr _1 , 
and fw is a correction factor which includes the most salient as- 
trophysical uncertainties. It is likely in the range 1—20, and we 



adopt a value of 10, following the arguments given in Cattaneo 
|etaL| ( |2009) . 

Typically we find kine tic luminosities of 5 — 10 x 10 45 erg 



with the calibration of Birzan et al. 



(2008), and -5 - 10 x 



10 46 erg s ! , corresponding to of order 0.01 — O.lJ^Edd with that 



of|Willott et al. ( 1999| l. Results for individual galaxies are listed 
in Table |2] 



5. Black-hole bulge mass relationship 

The scaling relationships of black holes and bulges are important 
tracers of the (possibly self-)regulated growth of black holes and 
galaxies across cosmic time. Accordingly, much observational 
effort has been invested in tracing the redshift evolution of these 
relationships (e.g . [Shields et al.||2003| [20061 1 Woo et al.||2"0"06] 
IPeng et al. 2006; Treu et al.|2007l IShen et al.||2008), although 

I 2 II i I II i I - II V ° 

sometimes with contradictory results. Overall these observations 
advocate a relatively mild redshift evolution of a few tenths of 
dex out to moderate redshifts (z<<1). 

We now compare our HzRGs with the local relationship, re- 
lying on our black-hole mass estimates ({ 4.3 and the photo- 



metric mass estimates of Seymour et al. (2007 1. Seymour et al 



estimated masses or upper limits for 70 HzRGs with Spitzer 
rest-frame near-infrared observations, constraining the relative 
contribution of the stellar component and dust at different tem- 
peratures. Photometric estimates of stellar masses are relatively 
robust for a range of star-formation histories, but poorly con- 
strained ages and initial mass functions cause systematic uncer- 
tainties of factors of up to about 2-3 for plausible ages between 
few x 10 8 yrs and about 10 9 yrs. 5 of our galaxies are included in 
the sample of |Seymour etal1 ( [2007l ). For TXS 1113- 178 which is 
not pa rt of this sample, we use the average mass of the [Seymour] 
et al.| sample, 10 115 M . We prefer to use this average mass 
rather than the K-bands magnitude derived from our SINFONI 
cube, because the latter may be significantly contaminated by 
the AGN, whereas the Seymour et al.| observations are near the 
rest-frame H-band where the stellar contribution to the overall 
spectral energy distribution is maximal. Moreover, the HzRGs 
of Seymour et al. ([2007 ) span a narrow mass range overall, not 
much greater than the uncertainties in the mass estimates of in- 
dividual galaxies. In the following we also assume that these are 
spheroid-dominated systems, consistent with the morphological 
studies of van Breugel et al. ( | 1998 ) and Pentericci et al. ( |2001 1, 
and plausible for high stellar masses of few 10 u M Q . 

In Figure |4] we show the black-hole bulge mass relationship 
of our HzRGs and compare it with the relationship for nearby 
galaxies d erived by|Haring & Rix| ( [2004[ ). Arrows mark galaxies 
for which Seymour et al. (2007i give upper limits on the stel- 



suggests that these upper limits are at most 20%-50% larger than 
the actual masses. This is negligible for our purposes. The sole 
exception may be MRC 1138-262 where the distinction between 
AGN and stellar continuum is more difficult due to the particu- 
larly powerful AGN. This analysis is based on IRAC and MIPS 
24/im and 60/zm photometry for all targets and includes a more 
systematic analysis of the dust emission (De Breuck et al. in 
prep.) 

We focus in the following on the ensemble properties of our 
sample, but each individual galaxy falls within the scatter of the 
local relationship (Figure [4]). This illustrates that HzRGs have at 
most about 4x larger black-hole masses for a given bulge mass 
than expected from the local relationship (a factor 2 if we use 
inclination-corrected black-hole mass estimates, see j ]4.3| ). This 
is broadly consistent with the redshift evolution of radio-loud 
AGN found by |McLure et al.| (|2006|) , if we extrapolate to the 
redshifts of our source s. lMcLure et al.l used the observed K-band 
magnitudes of radio galaxies to estimate bulge masses and steep- 
spectrum quasars out to z=1.5— 2.0, and predict an offset from 
the local relationship by 0.6 dex for z>2. 

The offset we find is significantly less than those reported 
for other populations of massive, strongly star-forming galaxies 
at z^2, (e g.|Walter et al.||2004[ [Alexander et al.||2008[ [CoppiiT 



|et al.|2008[|Wang et al.|2010[ ). The reason could be that we used 
stellar masses instead of dynamical masses based on the width 
of CO line emission like the above analyses, or it could indicate 
evolutionary effects. Redshifts z>l-2 are the most rapid phase 
in the evolution of massive galaxies, which makes observations 
at these redshifts particularly interesting, in spite of the observa- 
tional challenges. 

The offset of HzRGs relative to the local sequence is for- 
mally consistent with observations of the redshift evolution of 
the M-c relationship out to moderate (z~0.5) redshifts. Using, 
e.g., the relationship of | Woo et al. (2008 1, we would expect an 
offset by a factor ~ 2. However, it is not clear whether this is 
indeed an evolutionary effect or caused by selection effects. For 
example, our HzRGs, which are somewhat more massive than 
the local sample o f |Haring & Rix[ would fall on the relationship 
by Haring & Rix| ( |2004| ), if the latter was somewhat steepened 
(within their scatter). Alternatively, e.g., Wyifhe ( 2006[ ) propose 
that black hole and bulge properties could be related through a 
log-quadratic rather than a log-linear relationship, in which case 
our targets would show no sign of evolution. 

Observational biases may also contribute. For example, fol- 
lowing La uer et al.| ( 2007[ ) we may expect that the intrinsic scat- 
ter in the relationship combined with the the steep slope of the 
high-mass end of the galaxy mass function may bias the rela- 
tionship for massive galaxies towards higher black hole masses. 
We argue in ^7] that, whether this offset is intrinsic or not, mild 
offsets of factors 2-4 are not inconsistent with the basic picture 
where HzRGs are massive galaxies that have completed most of 
their growth by z^2. 



6. Low accretion rates of SMBHs in HzRGs 



In I 4.4 we used our black-hole mass estimates and broad-line 



lar mass. A new, ongoing analysis of the |Seymour et al.] sample 



luminosities to estimate accretion luminosities for our targets in 
the range of a few and up to about 10% Eddington. When reach- 
ing similar accretion luminosities of a few percent Eddington, 
Galactic "miniquasars" undergo rapid changes in their accretion 
properties, with dramatic effects for their X-ray and radio emis- 
sion (e.g., Nowak||1995 Maccarone et al.||2~003[ |Kording et al. 
2006a b i Fender et al. 2004J. The analogy with these stellar-mass 
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black holes has very recently inspired several authors to pro- 
pose a fairly complex scenario of "synthesis" ( |Merloni & H einz 
|2008] i or "unified" pmidakis et al.|2009] > AGN feedback, where 
an initial, strongly radiative phase with high accretion rates is 
followed by a later phase with low accretion rates, during which 
the radio source inhibits gas cooling and accretion through the 
injection of kinetic energy into the ambient medium (see also, 
e.g., Sijacki et al.||2007 1. Churazov et al.| (2Q05) proposes that 
this may be accompanied by episodes of highly unstable accre- 
tion, rapid "radio flaring", and powerful radio activity. |Merloni 
|& Heinz| ( |2008| ) use the local black-hole mass function and X- 
ray constraints on the radiative energy output of AGN to quantify 
the accretion and feedback history of AGN within this scenario 
in a semi-empirical way. 

The HzRGs of our sample have luminous > 10 46 erg 



nosities >0.01 J^Edd, similar to those in HzRGs, making this 



s l ) AGN as well a s very short-lived (a few x 10 6 7 yrs;|Blundell 



& Rawlings||1999[), powerful radio sources of up to 10 W 



Hz J at 500 MHz ( |Miley & De Breuck|[2008| , and they show 



the signatures of fast, AGN-driven kpc-scale outflows of sig- 
nificant fractions of their interstellar medium (Nesv adba et al.| 



2006 2008b ). These are several of the characteristics expected 



for the 'flaring' phase between radiative and kinetic feedback 
in the 'synthesis' models, and the low accretion rates of their 
black holes could be another. Churazov et al. (2005) postulate 
that during this phase, kinetic feedback and gas heating become 
more efficient as accretion rates drop further, until AGN and host 
galaxy reach a state where very low accretion rates can effi- 
ciently balance gas cooling ( |Croton et al.|2006] ). This may be the 
case in nearby, massive, early-type galaxies with radio sources 
(Best et al. 2006; Hardc astle et al.|2007| ) where star formation is 
significantly less efficient than in ordinary star-forming galaxies 
( [Nesvadba et al.|2010| ). 

For a more quantitative comparison we show the ratio of 



kinetic to radiative energy output (derived in {A.5 and S4.4 
respectively) in Figure [5] as a function of accretion luminos- 
ity. We adopted this Figure from Merloni & Heinz ( 2008 , their 



Figure 1.), and use the kinetic luminosities derived with the rela- 
tionship of |Willott etaLl ( |1999| ) and |Birzan et al.| ( |2008| ),respec- 
tively. For accretion luminosities above the critical value, the line 
splits to accommodate radio-quiet (lower branch) and radio-loud 
(upper branch) AGN. All of our targets are in the critical accre- 
tion regime, and have ratios of kinetic-to-radiative power that are 
near the values suggested by Merloni & Heinz (2008 ). 

Figure [5] also shows other populations of powerful AGN 
in massive, intensely star-forming galaxies at z>2, like sub- 
millimeter galaxies with BLR measurements ( [Alexander et al.| 
2008), and quasars at redshift z~2 and z~6, respectively (taken 



from |Coppin et al. 2008 ; Walt er et al.|2004||Wang et al.|2010| ). To 
estimate the kinetic luminosity of the radio sources in these sam- 
ples, we used the approach outlined in f 4.5 and published radio 
fluxes, where we assumed a radio spectral index of a — —0.8. 
For the high-redshift AGN, this has the character of an up- 
per limit, because the observed faint radio fluxes may be sig- 
nificantly contaminated, or even dominated, by star formation. 
Unlike the HzRGs, these galaxies do not fall into the critical ac- 
cretion regime. 

Similarly, the black holes in our HzRGs seem to have lower 
accretion luminosities than equally massive black holes studied 
in large quasar surveys at z^2-3. For example, Kollmeier et al. 
(2006 ) find that typical accretion luminosities in their X-ray and 
infrared-selected quasar sample taken from the AGES survey are 
a factor 10 higher (but tails of their distribution extend to the 
same range found for HzRGs). Their selection suggests >70% 
completeness for black holes with masses and accretion lumi- 



comparison robust (but see, e.g. Kelly et al. 2010, Schulze & 
|Wisotzki|2010| for a discussion of the effects of incompleteness 
in statistical studies of the accretion history of AGN). 

7. Implications for high-redshift galaxy evolution 

Observational evidence is growing that the host galaxies of 
HzRGs may be massive galaxies seen during a late stage of 
their active formation phase: At all redshifts, radio galaxies are 
among the brightest galaxies in the K-band, suggesting they are 
at the upper end of the mass function (De Breuck et al. [2001 



Willo tt et al."|200"3l see also |Seymour et al.||2007| ). At redshifts 



z> 3, most HzRGs are luminous submillimeter emitters with in 
tense star formation, whereas only few galaxies are detected at 
lower redshifts ( |Archibald efaLl[200T] |Reuland et aH|2004l ). If 
this is not entirely from selection effects, then it may suggest 
that HzRGs as an ensemble are transiting from actively star- 
forming to passively evolving during this epoch, perhaps as a 
consequence of AGN feedback signaled by the extended out- 
flows of warm ionized gas. It is interesting that this happens dur- 
ing the epoch where the red sequence seems to emerge ( Kodama 
|et al.|2007D . 

The results of §S]5] and [6] suggest that the black holes in 
HzRGs may also be nearing the end of their period of active 
growth. Compared to submillimeter galaxies and quasars at z>2 
([Alexander et al.||2008[ [Coppin et al.|[2008l | Walter et al.||2004 
Wang et al. 2010), they fall remarkably close to the local rela- 
tionship between black hole and bulge mass, as can be expected 
if the evolution of black holes and bulges is driven by the same 
mechanism. 

Figure|5]suggests that the accretion rates of HzRGs may have 
dropped to a critical value where mechanical feedback begins to 
dominate radiative feedback, unlike the black holes in the above 
mentioned submillimeter galaxies and quasars, which all have 
higher accretion rates. It is tempting to speculate that this tran- 
sition may also have triggered the large-scale outflows of the 
HzRGs, which may sweep up and heat the ambient cold gas out 
to scales of 10s of kpc. Integral-field spectroscopy of luminous, 
infrared-selected quasars with moderate radio sources reveals 
disturbed gas kinematics, but no evidence of equally strong out- 
flows extending over scales of 10s of kpc (Alexander et al. 2010 
Polletta et al., in prep.). 

Obviously, radiatively efficient accretion is maintained in 
HzRGs, with J% / =fewxl0 46 erg s , but outflow rates of a 
fewlOOMgyr-'of warm ionized gas (ignoring other gas phases 
Nesvadba et al. 2006 , 2008b ) suggest that the gas supply for star 
formation and AGN activity is being depleted rapidly. Broad UV 
absorption lines in MRC2025-218 ( |Humphrey et al.|2008| ), remi- 
niscent of broad absorption lines sometimes observed in quasars, 
suggest that outflows are also present near the black hole, mak- 
ing it plausible that, after the radio-loud phase, neither the su- 
permassive black holes in HzRGs nor their host galaxies, will 
have enough "fuel" to maintain much more growth, unless sig- 
nificant amounts of cold gas are supplied, for example through a 
major merger. The environments of HzRGs are rich in ionized, 
neutral, an d molecular gas ([V illar-Ma rtm et al.||2003 van Ojik 
|et al.|1997} |Nesvadba et al.|2009| ), but the analogy with nearby 
massive galaxy clusters illustrates that much of this gas can be 
heated by the AGN (see, e.g., |McNamara & Nulsen]|2007| and 
references therein). 

Moderate, subsequent bulge growth may occur through ac- 
cretion of satellite galaxies of HzRGs ( Nesvadba et al.| |2008b 



|Hatch et~a l. 2009). Ram-pressure stripping, tidal effects, and 



Nesvadba et al.: The black holes of HzRGs 



7 



on-going star formation will likely have removed much of the 
gas in the satellites by the time they merge, in which case 
galaxy growth would not necessarily lead to black-hole growth 
(Kazantzi dis et al.| |2005 ). This is consistent with structural pa- 
rameters of nearby massive early-type and cD galaxies ( jWeinzirl 
|et al.| |2009; Schombert 1987 ) and could eliminate the small off- 
set of HzRGs from the local black-hole bulge mass relationship 
(if it is not entirely due to systematic effects), or else change the 
position of our galaxies by a few dex within the large scatter of 
this relationship. 



8. Summary 

We presented an analysis of the black-hole properties in 6 pow- 
erful radio galaxies at z^2 (HzRGs), which show broad, lumi- 
nous Ha line emission that we identify as emission from the 
AGN broad-line region (HaBLR). This study is part of an on- 
going, systematic analysis of 50 HzRGs with rest-frame optical 
integral-field spectroscopy and with the goal of quantifying the 
role of AGN feedback for galaxy evolution at high redshift. Our 
main results are as follows. 

(1) HaBLRs in HzRGs are not rare. 6/31 HzRGs at redshifts 
where Ha can be observed show HaBLRs, corresponding to 
^20% of our sample. The broad Ha lines have emission-line 
luminosities of a few x 10 44 erg s _1 and line widths of FWHM> 
10,000 km s . Broad H/3 lines in 2 cases and X-ray constraints 
suggest low extinction (Ay ~ 1 mag). 

(2) We estimate black-hole masses and bolometric luminosi- 
ties from the Ha luminosities and line widths, finding masses 
of a few xlO 9 M and up to 10 10 M which is broadly con- 
sistent with the masses of the most massive black holes in the 
nearby universe. Bolometric luminosities are a few percent of 
the Eddington luminosity, lower than in other populations of 
massive, intensely star-forming galaxies and quasars with black 
holes of similar mass. 

(3) Black-hole and stellar masses suggest that HzRGs fall 
on the black-hole bulge mass relationship for local galaxies, for- 
mally with a small offset to lower bulge masses by at most a 
factor 4. Much of this offset may be from systematic effects; 
alternatively, it could be explained with an extended phase of 
accretion of (gas-poor) satellite galaxies. 

(4) Relative to the Eddington luminosities, HzRGs have ra- 
diative and (through the radio source) mechanical energy out- 
put rates that are in the "critical" range expected by recent "uni- 
fie'" AGN feedback scenarios. This could imply that they have 
reached the end of their period of active growth, if these models 
are fundamentally correct. We argue that the same is probably 
true for the host galaxies of HzRGs, suggesting that black hole 
and host galaxy could be terminating their phase of active growth 
at the same time. 
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Fig. 1. Spectra of 5 o f the HzRGs with BLRs seen in Ha and SINFONI observations (for the spectrum of MRC 1138-262 see 
Nesvadba et alj|2006] l. Top to bottom: MRC0156-252, MRC 10 17-220, TXS1113-178, MRC1558-003, and MR C2025-218. The 
BLR of MRC0 156-252 and MRC1558-003 were detected through longslit spectroscopy by |Larkin et al.| ( |2000| l and |Humphrey et al. 
(2007), respectively. All spectra are smoothed by 5 pixels (corresponding to 25A or ^340 km s l ) to emphasize the BLR. This is 
only for the sake of visualization, the analysis was performed on the unsmoothed spectra. 
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Fig. 2. Close-up of the Ha spectral region of MRC0156-003. We show the spectrum with the best fit superimposed (thick black 
line). Thin black lines show each individual fit component. The lower panel shows the fit residual. The line at A=2.03/im is the 
[SII]AA67 16,6731 doublet, where the two components are blended. At wavelengths shortward of 1.9jUm the atmospheric transmis- 
sion drops to <10%. 
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Fig. 3. We detected broad Hj3 line emission in two of our galaxies, MRC 1558-003 and MRC2025-218 (upper panels). Thick black 
lines show a Gaussian line profile with the same redshift and line width, which we measured for Ha, scaled to match the H/3 flux. 

( |2002| 



Thick gray lines show the Fell template from Iwamuro et al. 



to illustrate that the contribution of Fell is negligible in both 
cases. We subtracted the continuum in both spectra. The middle and lower panels show the fit residual spectrum after subtracting 
Hj3 and the position of night sky lines, respectively. 
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Fig. 4. Black-hole masses as estimated in f 4.3 as a function of stellar mass of the host galaxy, taken from Seymour et al.| (2007 i. 
Red dots mark the HzRG in which we detected HaBLRs, not correcting for inclination, and blue empty squares show the same 
galaxies with black-hole masses estimated assuming the inclination correction o f |Jarvis & McLu re (2002 ). The empty circle marks 
TXS 11 13-178 for which we used the average HzRG mass of Seymour et al. ( 2007| l. Arrows mark galaxies for which Seymour 
|et al.| only provide upper limits (see text for details). Empty circles show the nearby galaxies of |Haring & Rix| p004), and the black 
dashed line shows the best linear fit to their sample. 
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Fig. 5. The ratio of kinetic and radiative luminosity of the AGN as a function of accretion power, following 
(2008). Above a critical accretion power of a few percent of Eddington, the relationship splits into two branches to accommodate 
radio-loud and radio-quiet, rapidly accreting AGN with radiatively efficient accretion. Red filled and empty triangles show the radio 
galaxies of our sample, using different calibrations to estimate the kinetic energy from the observed radio power (§|4.5|[6| . Small 
blue squares and green circles show z^2 quasars and submillimeter galaxies from Coppin et al.| f2008 ) and Alexand er et al.| p008), 
respectively. Small light blue triangles show z>6 quasars taken from [Wang et al.| ( |2010| l. Blue arrows mark the typical accretion 
regime of massive early-type (and typically radio-active) elliptical galaxies. Only the HzRGs fall into the critical accretion regime, 



where (if the scenario of Me rloni & Heinz|2 008 and others is fundamentally correct) we may expect to find evolved supermassive 
black holes near the end of their active phase of growth, as they are transiting into the low-accretion regime spanned by massive 
early-type galaxies. 
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Table 1. Observed wavelengths, redshifts, full-width a half maximum, fluxes of the Ha and Hj3 broad lines, and V-band extinctions. 



Source 




Redshift 


FWHMfflaFl 


F(Haf 


1 


FfHfifl 


v \_\ 




[A] 




[kms- 1 ] 


[erg s~ 4 


J 


[erg s _T ] 


[mag] 


MRC0156-252 


19764.1 


2.01143 


12436 


1.75 x 10 


-14 






MRC1017-220 


18202.9 


1.77356 


12006 


5.3 x 10" 


15 






TXS1 113-178 


21253.7 


2.23851 


11063 


1.4 x 10" 


14 






MRC1 138-262 


20817.9 


2.17201 


14900 


1.9 x 10" 


14 






MRC1558-003 


23002.5 


2.50488 


12425 


8.3 x 10" 


15 


2.1 x 10~ 15 


0.9 


MRC2025-218 


23829.9 


2.63094 


8023.6 


5.5 x 10- 


15 


1.3 x 1Q- 15 


1.0 



" corrected for instrumental resolution 

for the broad line 
for the broad line 

V-band extinction of the broad line region. 



Table 2. Broad Ha luminosities, black holes masses, Eddington and bolometric luminosities, jet kinetic luminosities, and stellar 
masses. 



Source 


j^(Hafl 


M B #] 




^bo0 










[erg 


i ~ 


[m T 


[erg 




[erg s^] 


[erg 


s-T 


[erg s _1 f 


[Ms] 


MRC0 156-252 


5.0 x 


10 44 


1 x 10 1U 


1 X 


1Q 48 


5 x 10 46 


5 x 


10 4i 


7xl0 46 


<1.6xl0 12 


MRC 1017-220 


1.2 x 


10 44 


5x 10 9 


6x 


10 47 


1 x 10 46 


4x 


10 45 


2xl0 46 


<5xl0 u 


TXS1 113-178 


5.3 x 


10 44 


9x 10 9 


1 X 


10 48 


5 x 10 46 


6x 


10 45 


7xl0 46 


N/A 


MRC 1138-262 


6.6 x 


10 44 


2x 10 10 


3x 


10 48 


6 x 10 46 


2x 


10 46 


3xl0 47 


<1.2xl0 12 


MRC1558-003 


4.2 x 


10 44 


1 x 10 10 


1 X 


10 48 


4 x 10 46 


9x 


10 45 


2xl0 47 


5xl0 n 


MRC2025-218 


3.1 x 


10 44 


3x 10 9 


4x 


10 47 


3 x 10 46 


1 X 


10 46 


lxlO 47 


4xlO u 



a not corrected for extinction 

b masses corrected for a putative orientation bias would be factors 2 lower 
Eddington luminosities corrected for a putative orientation bias would be factors 2 lower, 
not corrected for extinction 



estimated with the relationship of BTrzan et al. 



(2008 



estimated with the relationship of |Willott et al. 
8 from Seymour et al. ( 2007 



1999) 



